Supercooled Liquid Dynamics Studied via Shear-Mechanical Spectroscopy 
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We report dynamical shear-modulus measurements for five glass-forming liquids (pentaphenyl 
trimethyl trisiloxane, diethyl phthalate, dibutyl phthalate, 1,2- propanediol, and m-touluidine). The 
shear-mechanical spectra are obtained by the piezoelectric shear-modulus gauge (PSG) method. 
This technique allows one to measure the shear modulus (10 5 — 10 10 Pa) of the liquid within a 
frequency range from 1 mHz to 10 kHz. We analyze the frequency- dependent response functions 
to investigate whether time-temperature superposition (TTS) is obeyed. We also study the shear- 
modulus loss-peak position and its high-frequency part. It has been suggested that when TTS 
applies, the high-frequency side of the imaginary part of the dielectric response decreases like a 
power law of the frequency with an exponent —1/2. This conjecture is analyzed on the basis 
of the shear mechanical data. We find that TTS is obeyed for pentaphenyl trimethyl trisiloxane 
and in 1,2-propanediol while in the remaining liquids evidence of a mechanical (5 process is found. 
Although the the high-frequency power law behavior uj~ a of the shear-loss may approach a limiting 
value of a = 0.5 when lowering the temperature, we find that the exponent lies systematically 
above this value (around 0.4). For the two liquids without beta relaxation (pentaphenyl trimethyl 
trisiloxane and 1,2-propanediol) we also test the shoving model prediction, according to which the 
the relaxation-time activation energy is proportional to the instantaneous shear modulus. We find 
that the data are well described by this model. 



I. INTRODUCTION 

The nature of the relaxation processes taking place in 
supercooled liquids approaching the glass transition has 
been a major subject of study for a number of years. 
Understanding how the different response functions are 
connected in such systems is still a fundamental goal to 
reach. It is not clear whether the various observables dis- 
play some universal features approaching the glass tran- 
sition. Although dielectric spectroscopy is the most com- 
mon experimental tool, a more detailed characterization 
of their behavior can be obtained by measuring differ- 
ent quantities, like the shear modulus, that are impor- 
tant from a practical as well as theoretical point of view. 
Although this dynamic variable can be measured above 
~ 10 2 Hz with conventional techniques, high-frequency 
data are scarce in the literature. 

Motivated by these reasons we employed the 
piezoelectric-shear-modulus-gauge (PSG) [l| method to 
measure the shear modulus of five glass-forming liquids. 
This technique allows us to measure the dynamic shear 
modulus of the supercooled liquids just above the glass 
transition where it takes values between 0.1 MPa and 10 
GPa. By means of the PSG technique we can easily ob- 
serve the a relaxation process in the shear response. The 
frequency range of the technique is wide (10~ 3 — 10 4 Hz), 
and we also observe a mechanical Johari-Goldstain (3 re- 
laxation 0] • As shown before 0] indeed, this technique is 
sensitive to the secondary process and we find evidence 
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of the presence of a shear (3 relaxation 0, 0, [f| in some 
of the mechanical spectra reported here. 

In Section II we describe the experiment performed 
and the liquids studied reporting the frequency resolved 
mechanical spectra. In Sec. Ill we report the main find- 
ings of our study. In Section III we present the analyzed 
data showing the temperature dependence of the shear- 
mechanical a relaxation frequency and a test of some 
conjectures and models about the dynamics. Finally we 
draw some general conclusions in Sec. IV. 



II. EXPERIMENT 

The piezoelectric shear-modulus gauge (PSG) method 
is based on the piezoelectric properties of the material 
that composes the measuring device. The piezoelectric 
transducer is formed by three discs made of a special ce- 
ramic compound jy] that has a pronounced piezoelectric 
effect (see Fig. [T]). The working principle of the PSG 
is illustrated in the inset of Fig. Q] where the one disc 
equivalent of the three discs system is shown. The ce- 
ramic disc is covered with a silver layer on both faces 
constituting the electrodes. When a voltage is applied, 
each disc expands or contracts depending on its intrinsic 
polarization and on the direction of the acting electric 
field. The electric capacitance of the disc depends on its 
strain state so that, if a material is partially clamping 
its motion, the measured capacitance will be lower than 
that of the free moving disc. Measuring accurately this 
electrical capacitance [7( we can obtain the stiffness of 
the medium adhering to the disc. In other words we can 
convert the electric impedance into the shear modulus 
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FIG. 1: (Main panel) Pictorial representation of the shear- 
transducer (PSG) during its expansion-contraction. The liq- 
uid sandwiched between the discs clamps their movement in- 
ducing a measurable change in the piezoelectric capacitance 
of the disc-system. The shear modulus can be deduced from 
this electrical capacitance. (Inset) The one-disc equivalent 
system of the PSG. The applied electric field causes an ex- 
pansion or a contraction depending on thepolarization of the 
disc (represented by the arrow). See Ref. [lj for more details. 



knowing the exact relationship between the two [l] . The 
three-discs geometry, used in the experiment (main panel 
of Fig. [lj is employed to reduce unwanted effect like the 
bending motion that would be present in a one-disc de- 
vice. The interested reader can find details about the 
technique in Ref. 

The measurements are performed cooling the liquids 
via a home-built closed-cycle cryostat Q. This has an 
absolute uncertainty on the temperature that is less than 
0.2 K and a temperature stability better than 2 mK. 

The liquids studied are the following: pentaphenyl 
trimethyl trisiloxane (DC705), dibutyl phthalate (DBP), 
diethyl phthalate (DEP), 1,2-propanediol, and m- 
touluidinc. The DC705 is a diffusion-pump oil from Dow 
Corning, all the other liquids were acquired from Sigma- 
Aldrich. No filtration or purification of the samples was 
performed before the measurement. 

All mechanical spectra are reported in Figure O Here 
the real part G' and the imaginary part G" of the com- 
plex shear modulus G{v,T) are presented^. The reac- 
tive and the absorbitive part of the shear-response of the 
liquids studied are measured at several temperatures. 

For the shear loss peak of the DC705 there is a clearly 
defined shear a process that shifts to lower frequencies as 
the temperature decreases. We have found no sign of a 
shear (3 relaxation in this liquid. This is also confirmed by 
the fact that the curve maintains the same shape lowering 
the temperature as discussed below. We will see also that 
the high-frequency side of the shear-mechanical spectrum 
shows a pretty constant slope. 



The shear response of DBP clearly shows a (3 relax- 
ation. The shape of the imaginary part of the response 
function is strongly deformed at high temperatures where 
the a process is merged with the (3. When the temper- 
ature is low enough, the high-frequency tail of the re- 
sponse functions shows a pronounced increase signaling 
an emerging [3 relaxation. Unfortunately the f3 loss peak 
lies at much higher frequencies than the upper limit of 
our device. The dielectric (3 process of this liquid has 
been observed [rJ to be around 10 6 Hz in a temper- 
ature window similar to the one of our mechanical ex- 
periment. This is consistent with the fact that the f3 
relaxation is outside the PSG dynamical range since the 
shear-loss peak frequency is generally higher than the 
dielectric one 0], (usually the two peaks are within the 
same decade). For future studies on this liquid the upper 
frequency limit of the PSG could possibly be increased 
by slightly modifying the transducer geometry. 

Evidence of a secondary loss peak is also found in the 
mechanical response of the DEP. The dielectric (3 relax- 
ation process in DEP has been the subject of an accurate 
investigation in connection with the behavior of the en- 
tropy in the supercooled and in the glass state [ll| , rais- 
ing fundamental questions about the nature of fast molec- 
ular motions in ultra-viscous liquids. Very recently this 
dielectric [3 process has been shown to be intimately re- 
lated to the results of positron-annihilation lifetime spec- 
troscopy We note that in those dielectric measure- 
ments the (3 loss peak is found at about 2 x 10 4 Hz at 
T = 179 K consistent with the fact that mechanically this 
process is above the upper limit of the PSG. We shall see 
in Sec. Ill that when the shear loss is plotted on a log- 
arithmic scale, the low-frequency tail of the (3 relaxation 
is more clearly visible. 

The spectra taken for 1,2-propanediol are character- 
ized by an high shear-modulus and no signature of (3 
relaxation is found. This hydrogen-bonded liquid has 
recently been studied in dielectric experiments and com- 
pared to polymers differing in chain length and in the 
number of OH group [l3T |. These dielectric experiments 
did not reveal a visible (3 process, and a temperature- 
independent Cole-Davidson stretching parameter was 
found. These findings are fully in agreement with our 
shear spectra since the shape of the response function is 
very weakly temperature dependent as it will be stressed 
in the following. 

Rather interesting mechanical spectra are found for m- 
toluidine. A weak, but non-zero high-frequency tail of the 
loss response can be noted in Fig. [2] This may hint to the 
existence of a shear (3 relaxation at frequencies around 1 
MHz, consistent with recent broad-band dielectric exper- 
iments performed on this liquid 1J] where a (3 relaxation 
process could be resolved, but only in the glassy regime 
where the relaxation frequency reaches a value around 1 
kHz at about 130 K. This could be a range in which this 
(3 process could start to be detected mechanically 15 1 . 
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FIG. 2: Real (left) and imaginary (right) part of the shear response of the liquids studied. (Form top to bottom) spectra of 
DC705 taken at the temperatures (from right to left): 249, 246, 243, 240, 237, 234, 231 K; of DBP at 188, 186, 184, 182, 180, 
178, 176 K; of DEP at 199, 197, 195, 193, 191, 189, 187, 185, 183, 181 K; of 1,2-Propanediol at 195, 192, 189, 186, 183, 180, 
177, 173, 171 K and of m-toluidine at 198, 196, 194, 192, 190, 188, 186 K. 
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FIG. 3: Frequency of the mechanical a peak as a function 
of the temperature for DC705 (O, T g = 224.0 K), DBP ( □, 
T g = 172.3 K), DEP (*, T g = 179.5 K), 1,2-Propanediol (0, 
T g = 157.8 K), m-toluidine (A, T g = 181 K). 



Finally we mention that we attempted to mea- 
sure the frequency-resolved shear-modulus also for 1,3- 
propanediol and for propylene carbonate. An anomalous 
shear response was found for these two liquids, however, 
signaling a probable crystallization. 



III. DISCUSSION 

The first information that can be extracted from our 
measurement is the a process loss peak frequency v max . 
This has been deduced from the imaginary part through 
a simple unbiased method. We fitted the closest six 
points to the maximum of log 10 (G"/ GPa) as a function 
of lo 

§io (^/Hz) with a second order polynomial to iden- 
tify the maximum G max and the corresponding frequency 
Vmax ■ These loss-peak frequencies are reported in Fig. [3] 
as function of the scaled temperature (T — T g ). Here we 
define T g as the temperature where the v max = 1 mHz. 
T g is identified by a linear extrapolation of the last three 
points of log 10 (V/Hz) as a function of T. 

Time-temperature superposition (TTS) is a property 
of the a process appearing in the susceptibilities of some 
super-cooled liquids expressing the fact that the shape of 
the response function remains the same when the system 
is cooled. When TTS applies, the response is only shift- 
ing its characteristic time and its amplitude if we cool the 
liquid. Mathematically speaking the generic response x 
can be rewritten as T) = A(T)$(v / v max (T)) where 
(T) is a characteristic frequency depending on T. 
Having found the two parameters G max and v max at ev- 
ery given temperature, it is rather easy to check wether 
TTS applies for a given liquid. 

In Fig. 0] we report the imaginary part of the shear 
modulus divided by G max as function of vjv max . This 



plot underlines that for some of the liquids (DC705 and 
1,2-propanediol) the shape of the alpha process remains 
constant lowering T. This is not the case for DBP and 
DEP where the curves are detaching at high frequencies. 
The lack of TTS is clearly related to the presence of a 
(5 peak, indeed the logarithmic scale used in this plots 
makes the the low- frequency tail of the secondary process 
more visible. For m-toluidine the situation is ambiguous. 
Indeed for this liquid TTS seems to apply, but if we look 
at the high frequency end on a smaller scale (Fig. [4] right- 
bottom) there is a sensitive difference between the curves. 
In this figure we report also the scaled spectra of DC705 
that satisfy TTS on the small scale for an easier compar- 
ison. We also note that the dielectric strength of the sec- 
ondary relaxation reported for m-toluidine [3] is rather 
low while in the shear-mechanical spectra it seems to be 
more intense in comparison with the a process. This 
is another important feature of our mechanical measure- 
ment: a weak dielectric secondary process is magnified in 
the shear response [H, 0| ■ 

An interesting feature of the spectra of 1,2-propanediol 
is underlined by the TTS plot in Fig. |4j The portion 
of the spectrum at frequencies lower than v max (i.e. at 
v/v-max < 1) is slightly deformed. The shape of this re- 
sponse function is reproducible and it seems to be almost 
temperature-independent in the temperature range ex- 
plored in our experiment. This "shoulder" could be con- 
sidered a peculiar feature of the mechanical relaxation of 
1,2-propanediol since it is not found in dielectric spectra 
(see for example Ref. Note that an additional pro- 

cess (with similar shape) close to the alpha relaxation is 
found in the case of tri-propylene glycol in high pressure 
conditions (see Ref. [l6(). Further studies of the dielec- 
tric relaxation of 1,2-propanediol at high pressure would 
be needed to better understand the origin of the shape 
of its mechanical response. 

We want to stress again that, from our study, the lack 
of TTS seems to be intrinsically related to the presence 
of a secondary mechanical process. All the liquids stud- 
ied here that show the (3 relaxation (as also confirmed 
by the dielectric measurements cited in Sec. II), display 
a clear temperature-dependence in the shape of the me- 
chanical response function. Anyway it is not generally 
clear how the presence of the beta relaxation is related 
to the nature of a specific liquid. The understanding of 
the validity of TTS in some supercooled liquids requires 
the understanding of one of the most puzzling phenomena 
of ultra- viscous liquids and glasses that is the secondary 
relaxation process. 

It has been suggested from empirical observations [13, 
[IH that when TTS applies, the high frequency decay of 
the dielectric loss is characterized by a power-law behav- 
ior v a with exponent a = —1/2. A theoretical expla- 
nation for this power-law dependence of the dielectric 
response has been proposed assuming the dominance of 
long- wavelength fluctuations TjJ 2^, 21 1 . We have tested 
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FIG. 4: Shear loss G" of liquids studied scaled by the maximum loss G^nx as a function of the frequency v divided by the 
frequency of the peak v max (see also Fig. DC705, DBP, DEP, 1,2-propanediol and m-toluidine. In the right-bottom figure 
we show a comparison between the scaled spectra of DC705 (full line, full symbols) and m-toluidine (dashed line, open symbols) 
on a smaller scale. 



this conjecture for the mechanical response finding the 
minimum value of the logarithmic derivative of G" , i.e. 
ctmin = \d\nG" Jdhxv\ m i n . It has to be noted that the 
extraction of this information is much more complicated 
for the shear loss than the dielectric response. This is due 
to the fact that the shear response is generally affected 
by a higher noise level than dielectric. An example of the 
procedure followed to extract a m in is reported in Fig. [5] 
for the mechanical response of m-toluidine at T — 186K. 
When a is plotted as a function of the frequency a min- 
imum can be identified within some data points. The 
average of these points is taken as the a m i n . Note that 
a Debye-process would have a characteristic behavior of 
the logarithmic slope showing a sudden drop of a at the 
relaxation frequency v max . 



In Fig. we report the value of a m i n as a function 
of the loss peak frequency for the temperatures where 
the minimum of the derivative of the logarithm was well 
defined. We note that even if the behavior of this quan- 
tity for the liquids studied is compatible with a limiting 
a = —1/2 at low-temperatures, the data seems to lie sys- 
tematically above the —0.5 line. We stress also that is 
easy to expect that an v~ x l 2 behavior in the shear loss is 
disturbed more when a (3 relaxation is present since its 
intensity is enhanced in the mechanical spectrum com- 
pared to the dielectric. 

Finally we remark that the parameter ot m in is not ex- 
pected to monotonically decrease in a wider temperature 
range. Indeed at higher temperatures the shape of the 
shear-response should recover a simpler Debye-like shape 
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FIG. 5: Representation of the procedure used to find the mini- 
mum logarithmic slope a m i n . In this figure we report the anal- 
ysis of the mechanical spectrum of m-toluidine at T = 186K. 
The derivative of log(G"/GPa) with respect to log(V/Hz) is 
plotted as a function the frequency (stars). Some points are 
selected (red circles) and their average is taken to calculate 
amin (red line). Note that a Debye-like mechanical process 
with the same relaxation time would have a sharper drop of 
the a parameter reaching —1 at high frequencies. 



-0.1 
-0.2 
s i -0.3 
-0.4 
-0.5 



on the assumption that the relaxation in the supercooled 
liquid takes place with a local volume increase and that 
the activation energy is mainly elastic energy spent to 
shove aside the surrounding of the rearranging molecules. 
As it has been reported this model can be derived esti- 
mating the barrier height in a classical energy landscape 
approach [2^|. 

The shoving model predicts that the relaxation time 
is related to the infinite-frequency shear modulus Goo by 
the equation: 



ln(r) = In r + V c G ao (T)/T 



(1) 



where V c is a characteristic volume in a relaxation process 
(V c is assumed temperature independent) and Tq 1 is the 
phonon- frequency (tq ~ 10~ 14 s). The determination of 
Goo is not necessary to test the model if TTS applies, 
indeed if this property holds we can write 



(2) 



since the constant factor A determining the amplitude 
of the complex response in x(i/, T) — A{T)<&(v / v max {T)) 
is the same for the real and the imaginary part <£' and 
Note that this proportionality does not apply when 
TTS does not hold. For example in a case where the /3 
relaxation is present the value of Goo is modified by the 
secondary process and a more complicated fitting pro- 
cedure would be needed to estimate the value of the 
infinite-frequency shear-modulus. On the other hand, 
when TTS is satisfied, we can write the simple equation 
(using Eq.© in Eq.Jlj)) 



iogioM =logi To 



max 



(T)/T 



(3) 



-3-2-10 12 3 4 

og,„(v /Hz 

a 10* max ' 

FIG. 6: Minimum value of the logarithmic slope (a — 
din G" /din v) for the liquids studied: DC705 (O), DBP ( 
□ ), DEP (*), 1,2-Propanediol «>), m-toluidine (A). 



where the stretching parameters (the KKW (3 parameter 
for example) come close to the unity [22]. The relaxation 
frequencies belonging to this high temperature regime 
are currently outside the dynamical window of the PSG 
technique [23|. 

To conclude the analysis of our data we want to present 
a comparison between the mechanical spectral features 
and a simple phenomenological model for the dynamics 
of supercooled liquids. The shoving model [13] is based 



where B is a constant factor. Eq (J3j) represents the pre- 
diction of the shoving model when TTS also applies: the 
logarithm of the relaxation time is a linear function of 
the quantity G'^ nax {T)/T. This is tested in Fig. for 
the liquids in which time-temperature superposition is 
found to hold (DC705 and 1,2-propanediol). In this fig- 
ure the relaxation time (defined as r = (27ri/„ laa; ) _1 ) 
is reported as a function of l/T and as a function of 
Gm ax /T To report all the data in the same plot the 
abscissa has been normalized to the unity as T g /T and 
as x = {T g /G'^ ax J(G^ ax /T). Note that here we find 
G m ax,g through a linear extrapolation of the last val- 
ues G^ aa , to T g previously identified to give v ma x(T g ) = 
10" 3 Hz (corresponding to r(T g ) = (2tt10 3 )~ 1 ~ 150s). 

The dashed line in Fig. [7] represents the shoving model 
prediction (no adjustable paremeters are used in this 
function) ending, at high temperature, at the physically 
reasonable prefactot 10~ 14 s. 
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FIG. 7: Characteristic time r = l/^ ma i of the shear loss as 
a function of the scaled temperature T/T g (full symbols) and 
as function of x — (T g /G'J naWi g)(T/G'J nax ) (open symbols) for 
the liquids in which TTS applies: DC705 (O, T g = 224.0 K) 
and 1,2- Propanediol (<[>, T g = 157.8 K). The dashed line is an 
exponential function with prefactor to = 10~ 14 s and taking 
the value r = (27rl0' j )~ 1 ~ 150s s for x=l (no free parameters 
are adjusted in this function). 



IV. CONCLUSION 

We have reported shear-mechanical spectra of five 
glass- forming liquids close to T g . Via the PSG [l[ tech- 
nique we have investigated the behavior of the mechan- 
ical a process and found evidence of the presence of a 
mechanical /3 relaxation in dibutyl phthalate, diethyl ph- 
thalate and in m-toluidine. Time-temperature superposi- 
tion for the mechanical susceptibility is found to hold for 
the liquids without signature of j3 process (pentaphenyl 
trimethyl trisiloxane and 1,2-propancdiol). The conjec- 
ture originally developed for the dielectric response that, 
when TTS applies, a v~ x l 2 decay is found for the high fre- 
quency part of the loss is checked for the shear response. 
We find that the data are consistent with a limiting u^ 1 / 2 
low-temperature behavior although the minimum slopes 
are systematically higher than —0.5. The shoving model 
has been tested for the two liquids without (3 relaxation 
finding that it well describes the experimental data. 
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